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In this work we propose a definition of an Euroattractor. an attracting invariant measure of a 
certain iterated functions system (IFS). An IFS is defined by specifying a set of functions, defined 
in subsets of or in a classical phase space, which act randomly on the initial point, so it may 
be considered as a generalization of the notion of classical dynamical system. If the functions are 
sufficiently contracting, there exists an invariant measure of the system, often concentrated on a 
' fractal set. We investigate invariant measures of a certain class of generalized or weakly contracting 

O ' IFS's. 
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^ , I. INTRODUCTION 



Any function / which maps a phase space Q into itself may be considered as a classical dynamical system. An 
arbitrary point x £ fl determines the trajectory uniquely, x n = f n (xo), where the natural index n plays the role of 
time. Dynamical system defined by / is deterministic. There exist different ways to introduce randomness into the 
system. One possibility relays on taking into account an additive noise, x n +i = f(x n ) +£n, where £„ are uncorrelated 
random variables drawn according to a prescribed probability distribution P(£). The deterministic case is obtained 
Ph ■ if P(£) tends to the Dirac delta distribution, P(£) = 

In this work we analyze another class of stochastic systems, in which the choice of a system is random at each step. 
I . Such systems are called Iterated Function System (IFS), and may be considered as a generalization of the notion of 
' a classical dynamical swstem. 
7— I . An IFS is defined p| by a set of k functions ft : f2 — > Q, which represent dynamical systems in the classical phase 
space n. The functions ft act randomly with given probabilities pi. They characterize the likelihood of choosing a 
particular map at each step of the time evolution of the system. In general, the probabilities may be place dependent. 

Analyzing a given iterated function system, one may ask, how an initial point xo G is transformed by the random 
process. One may also consider a more general question, how does a probability measure /i on Q change with respect 
. to the Frobenius-Perron operator M associated with the IFS. If the functions /, are strongly contracting, then there 
*"^> ' exists unique invariant measure /i* of M - see for instance QL|^|^,Q and references therein. 

Interestingly, for a large class of IFS's the invariant measure /U* is characterized by a fractal structure. If the 
invariant measure is attracting, such IFS's may be used to generate fractal sets in the space fi. In particular, there 
exist iterated function systems leading to popular fractal sets: Cantor set, Sierpinski carpet or Sierpinski gasket [lj. 
In fact by considering a slightly generalized class of IFS's it is possible to design a random system, the invariant 
measure of which possesses any desired properties. In this work we shall visualize the statement, by presenting an 
Euroattractor - a system, the attractor of which has the shape resembling the contours of Europe. 

This paper is organized as follows. In the next section we recall the definition and basic properties of the standard 
IFS. Section III is devoted to a weakly contracting IFS, which exhibit nontrivial invariant measures, although they 
do not fulfill the standard contractivity assumptions. In section IV we define a generalized class of the systems and 
present the definition of an Euroattractor. 

II. ITERATED FUNCTION SYSTEMS 

An iterated function system (IFS) is defined by a set of k functions (dynamical systems) /,-, i = 1, ...,k, defined 
on a compact space and the set of k positive probabilities pi, which sum to unity. In the IFS of the first kind the 
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probabilities are constant, while in the more general, IFS of the second kind |J,|5j, they are functions of the position, 
Pi = Pi(x). For any x € CI the condition J2i=iPi( x ) = 1 nas to be fulfilled. In general, IFS T is defined by a set 
T = {O, /i : -> Q,pi : ft -> [0,1]} with i = l,...,fc. 

Any IFS generates a Frobenius-Perron operator M, which describes the time evolution of the probability densities 
on Q. In the simplest case, in which the space fl is an interval in K and the maps fi are invertible, any initial 
probability density p(x) is mapped into p'(x) by the Frobenius-Perron operator @ 

k 

P '(x) = m[ P (x)] = ^(frHxMfrH*)) 

i=l 

If p(x) is a non-singular probability density, so is p'{x). However, in the limit n — > oo the image M™[p(x)] needs not 
to be a smooth probability density, so it is advantageous to work in the larger space of probability measures on fi. 
This space includes also also singular measures like the Dirac delta distributions and their convex combinations. 
If an IFS is sufficiently contracting, i.e. 

a) the functions /, satisfy the Lipschitz condition, d(fi(x), fi(y)) < Lid(x, y)) for any x, y € fi. Here d(x, y) denotes 
the distance between both points and the Lipschitz constants Li are assumed to be smaller than unity, 

b) the probabilities Pi(x) are continuous functions, strictly positive for any x € fi, the IFS is called hyperbolic and 
there exists a unique invariant probability measure /x» satisfying the equation Mp* = p* |6|. 

Interestingly, the conditions a) and b) are not necessary to assure the existence of a unique invariant probability 
measure. Some other less restrictive sufficient assumptions which regard also IFS of the second kind, with place 
dependent probabilities, were analyzed in [2I 13, I5L I8M91 ITol Hll ] . 

The invariant measures of an IFS display often fractal properties. For instance, consider an IFS of the first kind 
defined on an interval and consisting of fc = 2 affine transformations, 

T x = {V = [0, 1], fi(x) = x/3, f 2 (x) = x/3 + 2/3; Pl = p 2 = ±}. (2.2) 

Since both functions are continuous contractions with Lipschitz constants L\ = L 2 = 1/3 < 1, and the probabilities 
Pi are constant and nonzero this IFS is hyperbolic. Thus there exists a unique, attracting invariant measure p*. It is 
not difficult to show pj that p* is concentrated uniformly on the Cantor set of the fractal dimension D = In 2/ In 3. 

Although /i» is singular, it is possible to give a general prescription, how to compute integrals with respect to this 
invariant measure. This is due to the fact that the measure /x* is attractive, i.e., M n p converges weakly to p* if 
n — » 00. Therefore, in order to obtain the exact value of an integral J n h(x) dp* for any continuous function h(x) it is 
sufficient to find the limit of the sequence h n := J n h(x) d(M n p) for an arbitrary initial measure p. This method of 
computing integrals over the invariant measure p* is purely deterministic |l|, |5|, lid Il2| . Alternatively, one may employ 
a random iterated algorithm by generating by the IFS a random sequence Xj £ fi, j = 0, . . . , n — 1 which originates 
from an arbitrary initial point .tq. Due to the the ergodic theorem for IFS's 0,0,0] the mean value ^Ej=o h{ x j)] 
converges in the limit n — > 00 to the desired integral, J n h(x) dp*, with probability one. 

To show an example of an IFS of the second kind let us change probabilities in T\ and define 

^ 2 = {fi=[0,l], fi(x)=x/3, f 2 (x)=x/3 + 2/3; Pl (x) = x, p 2 (x) = 1 - x}. (2.3) 

The probability P2 vanish at the point x = 1, so this IFS is not hyperbolic. In spite of this fact there exists the unique 
invariant measure p* concentrated in a non-uniform way on the Cantor set [5| • In this case displays multifractal 
properties, since the generalized dimension of the measure depends on the Renyi parameter [14j. 

III. WEAKLY CONTRACTING IFS 

Contraction conditions a) and b) imply existence of a unique invariant measure for the IFS. However, these condi- 
tions are too strong, and there exist IFSs with an attracting invariant measure which do not obey them. If the space 
il C M m and the functions fi are linear, the Lipschitz condition a) is equivalent to the statement that the moduli of 
all elements of the Jacobi matrices Ji :— dfi(x\, x m )/d(x\, ...,x m ) are strictly smaller then one. 

This is not the case for the following IFS defined on the square 

f n=[0,l]x[0,l], fc = 2; Pl = P2 = l/2 } 
= { h(x,y) = (xy, y) >, (3.1) 

[ f 2 (x,y) = ((x-l)y + l, y + S) J 



dx 



(2.1) 
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since the Jacobi matrices 

Ji=(l i) and J 2= (jj X ~ 1 ^j (3.2) 

contain entries equal to unity. This reflects the fact that both functions /i and /2 are not contracting in the y— 
direction. Numerical experiments suggest that in spite of this fact for irrational values of 5 -C 1, the associated 
Frobenius-Pcrron operator M possess a unique invariant measure /i*. Such an IFS, not fulfilling the condition a) or 
b) will be called weakly contracting. 
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FIG. 1: Invariant density of the weakly contracting IFS Tz ft~i.lfl obtained for 5 — \/E/W00. The cross-section of its support at 
y = 1/3, denoted by a dashed line, forms the Cantor set. 

Each random trajectory of the IFS Tj, exhibits a slow diffusion in the y-direction, the speed of which is governed by 
the rotation parameter S. Due to periodic boundary conditions, (x,y) = (ac| mo di, y|modi), every trajectory is confined 
to the square fl. If the parameter S is chosen to be irrational, there are no periodic trajectories. If 5 is sufficiently 
small than the contraction in the a;-direction will dominate the diffusion in y, so the trajectory will be sticked to the 
support of the invariant measure fj,*. Figure 1 shows a typical trajectory consisting of 6000 points generated by this 
IFS. To present the asymptotic properties of the invariant measure the first 100 points were omitted. 

Note that if the parameter 5 is set to zero, the variable y would be fixed to the initial value yo- Such an IFS would 
have an invariant measure concentrated on the interval y — yo. Analyzing the contraction constant we see that for 
yo = 1/3 its support is equal to the Cantor set generated by the Cantor IFS 1)2.2(1 . In general, for any yo < 1/2 its 
dimension would be 

, x In 2 , N 

D(yo) = • (3.3) 

Inyo 

We conjecture that this formula gives also the fractal dimension of the horizontal cross-section at y = yo of the 
invariant set of the IFS (|3.1(l in the limit 5 — > 0. 

Another example of an invariant measure of a nonlinear, weakly contracting IFS is presented in Fig. 2. The system 
is defined by 



O = [0,1] x [0,1], k = 

f\ : \ h{x,y) = 

h(x,y) = 



2; pi=p 2 = 1/2 
(iicsin 2 (7ry), y) 
(i(x-l)sin 2 (7ry) + l, y + S) 



(3.4) 
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As in the previous example each random realization of the IFS generates a trajectory confined in the torus, and for 
small values of the shift parameter S the contraction in the horizontal direction dominates the vertical diffusion. For 
(5 = the IFS reduces to one-dimensional, and the fractal dimension of its invariant set reads 



D(y ) = 



In 2 



In 2 — 2 lnsin(7r/?/o) 



(3.5) 



This number is smaller than one for all initial values apart of yo — 1/2. Numerical results performed for S <C 1 suggest 
that the same formula describes the dimension of the cross-section at y = yo of the invariant set of the IFS (|3.4fl in the 
limit (5 — > 0. In other words, the support of the invariant measure of the IFS T± may be considered as a collection 
of Cantor sets, with a continuously varying dimension D(y ), placed horizontally, one parallel to another. 




FIG. 2: Invariant density of the weakly contracting IFS 13.40 obtained for S — \/5/i000. 



IV. REPELLERS AND GENERALIZED IFS 



In the original definition of the IFS, presented in section II, all functions fi map the space f2 into itself. One may, 
however, relax this condition, and consider generalized IFS constructed of functions which map some fragments of f2 
out of this space. Such functions are called repellers, since they do not conserve the probability, and during the time 
evolution the probability flux leaks out of the system. 

A simple 1-D repeller is obtained from the standard logistic map, f(x) — rx{\ — x), if the control parameter r 
exceeds the critical value, r c = 4 |18| . In such systems a randomly selected initial value leads to a trajectory, which 
escapes the system with probability one. The complementary set of initial points which generate infinite trajectories 
is of measure zero and often forms a fractal set. The same property is characteristic to a large class of chaotic 1-D 
maps with a gap - a subset of £1 map_ped out of fi. An analysis of the fractal dimension of the invariant set for the 
tent map with a gap is presented in 19] . 

In this chapter we will consider an IFS defined on the square, = [0, 1] x [0, 1], consisting of k repelling functions. 
Assume that each function fi restricted to a given set Xj C acts as an affine transformation defined by real matrices 
Ji of size 2 and two components translation vectors Ki, but the complementary subset f2 \ X, is mapped out of O, 

X ) = { ' h ( V ) + X<EXl (4.1) 

V J I out of fl for x i Xi 
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Such functions fi are repelling, so instead of analyzing single trajectories, generated by a random initial point, it 
is convenient to work with initially smooth probability measures. Let us take for the initial measure /j,q the uniform 
(Lebesgue) measure on f2. To iterate the repelling IFS we find the images of the measures covering each of the subsets 
Xi by the maps fi, i = 1, k and add all contributions with prescribed weights - probabilities pi. More formally, we 
use the the dimensional analogue of the Frobenius-Perron operator 12.1fl 

k 

p'(x,y)=M[p(x,y)]=Y,Pi(f l rl (x,y))p{fr 1 &y)) 

i=l 

where the contributions from the function fi have to be taken into account only if the analyzed point (x, y) belongs 
to fi(Xi). Then the preimage f~ 1 (x, y) exists, so the i-th term in the expression i|4.2[l is meaningful. 




dfi(x,y) 



d{x,y) 



(4.2) 



FIG. 3: Generalized IFS with 'repelling' functions fi which send the polygons Xi into X[, and the remaining set Q. \ Xi is 
mapped out of the system. 

As an example of such a generalized, repelling system we consider the IFS T§ consisting of k = 13 functions with 
equal weights, pi = 1/13. The detailed form of the affine functions fi as well as the definitions of the sets Xi are 
provided in the Appendix. Figure 4 shows time evolution of the initially uniform measure fiQ. The right hand side 
of the figure shows the measures n n — M n \po], while the left hand side shows the sets, where the measure fi n is 
concentrated. Initially uniform measure fj, converges fast to the asymptotic measure [i*, and already after n = 6 
iterations can hardly be distinguished from it. The shape of the support of /i* explains why the name Euroattractor 
has been concocted for this system. 



V. CLOSING REMARKS 



The purpose of this work was to demonstrate usefulness of the concept of Iterated Function System. This general- 
ization of the concept of a dynamical system incorporates deterministic and stochastic behavior and may be applied in 
various fields of physics, and in general, of nonlinear science. A possible extension of the concept of IFS for quantum 
mechanics was recently worked out in [l5j . 

The IFS's belong to a larger class of random systems studied in (l6llri|. From a mathematical point of view several 
questions concerning the rigorous theory of IFS remain open. For instance, the sufficient conditions, under which 
an IFS possesses a unique, attracting invariant measure are known, but the formulation of the necessary conditions 
remains still an exciting mathematical challenge 0, IToL I20L l2l| . In this work we have introduced a weakly contracting 
IFS's Tz and T^, which do not satisfy the standard contractivity conditions, but do posses attracting invariant 
measures. 

For any given set S in the plane one may try to find an IFS, such that its invariant measure /i* has support on S. 
This versatility of IFS's allows one to apply them for image synthesis, image processing and image encoding 0,0- 
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FIG. 4: Euroattractor: the time evolution of the support (left) and the measure (right) of the initially uniform measure /io in 
the square. Already after n = 6 iterations the measure cannot be distinguished from the invariant measure fi* . 



For instance, the attractor of a generalized IFS described in this work, has the shape resembling the contours of 
Europe p^. 

It is a pleasure to thank Radoslawa Bach and Wojciech Slomczyriski for constructive remarks and fruitful interac- 
tion. This paper was presented at the European Interdisciplinary School on Nonlinear Dynamics Euroattractor 2002 
organized in Warsaw in June 2002. We are grateful to Wlodzimierz Klonowski for an inspiring name of the conference 
as well as for the opportunity to present this work during that event. 



APPENDIX A: AN EUROATTRACTOR 



In this appendix we provide an explicit definition of an Euroattractor - a generalized, repelling IFS !F§, the invariant 
measure of which is shown in Fig. 4. The space is equal to an rectangle. To address each pixel in a simple way its 
size is taken to be 500 x 600. 

The IFS consists of k = 13 affine functions fi defined in Q4.1J1 . and all the weights are constant, Pi(x,y) = 1/13. 
The parameters defining the corners of the quadrangles X^, and the elements of the transformation matrices Jj and 
the translation vectors Hi are collected in Table 1. Observe that some entries of the Jacobi matrix J§ have modulus 
larger than one, so this IFS does not satisfy the contractivity condition. 
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Table 1. Parameters defining an Euroattractor - the IFS T 5 



Ji 




220 
350 



300 
500 



0.37 0.66 
-0.74 0.66 




520 
430 



340 
320 



420 
530 



240 
420 



-0.33 
0.80 



-0.73 
-0.59 



626 
355 



170 

350 



265 
480 



230 
320 



325 
450 



0.52 0.38 
0.01 0.68 



16 
218 



340 
300 



240 
420 



500 
400 



400 
520 



-0.34 
-0.61 



0.54 
-0.51 



192 
814 



380 
500 



500 
300 



250 
450 



370 
250 



-0.17 
0.96 



-0.95 
0.08 



812 
31 



240 
460 



340 
550 



260 
440 



360 
530 



1.44 
-0.23 



-1.05 
1.26 




200 
490 



50 
590 



150 
270 




370 



-0.31 
-0.98 



-0.47 
0.22 



510 
430 



460 
540 



305 
330 



245 
440 



-0.47 
0.39 



-0.032 
0.99 



230 
280 



230 
450 



460 
280 



460 
450 



0.26 
-0.60 



1 

0.61 



10 



440 
280 



440 
450 



245 
280 



245 
450 



-0.51 
0.87 



0.70 
0.41 
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250 
280 



310 
190 



220 
190 



270 
100 



-0.22 
0.33 



-0.37 
0.11 



439 
456 



12 



330 




200 
190 



410 
100 



280 
290 



-0.60 
0.47 



-0.51 
-0.38 



359 
187 



13 



420 
10 



420 
90 



330 
10 



330 
90 




243 
30 
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